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Abstract. The QCDSF/UKQCD collaboration has an ongoing program to calculate nucleon matrix ele- 
ments with two flavours of dynamical 0{a) improved Wilson fermions. Here we present recent results on 
the electromagnetic form factors, the quark momentum fraction (a;) and the first three moments of the 
nucleon's spin-averaged and spin-dependent generalised parton distributions, including preliminary results 
with pion masses as low as 320 MeV. 

PACS. 12.38.Gc Lattice QCD calculations - 13.40.Gp Electromagnetic form factors 



1 Introduction 

The ability of generalised parton distributions (GPDs) 
to describe both exclusive and inclusive processes has led 
to an enormous amount of interest in these functions both 
experimentally and theoretically. Not only do GPDs en- 
compass the ordinary electromagnetic form factors and 
parton distribution functions, but they also allow for the 
computation of the total quark contribution to the nucleon 
spin [2] as well as revealing important information on the 
transverse structure of the nucleon A full mapping 

of the parameter space spanned by GPDs is an extremely 
extensive task which needs support from non-perturbative 
techniques like lattice simulations. 

Substantial progress has already been made in com- 
puting the first three moments of unpolarised, polarised 
PCTT l and tensor ^ GPDs on the lattice. 

In this paper we present recent results from the QCDSF 
/UKQCD collaboration. In section[21wc investigate the 
dependence of the Dirac and Pauli electromagnetic form 
factors, while section |3 contains preliminary results for 
the average fraction of the nucleon's momentum carried 
by the quarks, {x). Finally, in section0]wc present results 
for the first three moments of the GPDs H and H . 

2 Electromagnetic form factors 

The study of the electromagnetic properties of hadrons 
provides important insights into the non-perturbative struc- 
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ture of QCD. The EM form factors reveal important infor- 
mation on the internal structure of hadrons including their 
size, charge distribution and magnetisation. Phenomeno- 
logical interest in these form factors has been revived by 
recent Jefferson Lab polarisation experiments [5| measur- 
ing the ratio of the proton electric to magnetic form fac- 
tors, ^^^'^ G^i\q^) / G'>m {q^) ■ These experiments show that 
this ratio unexpectedly decreases almost linearly with in- 
creasing (7^, indicating that the proton's electric form fac- 
tor falls off faster than the magnetic form factor. 

A lattice calculation of the q^ dependence of the pro- 
ton's electromagnetic form factors can not only allow for a 
comparison with experiment, but also help in the under- 
standing of the asymptotic behaviour of these form fac- 
tors. Such a lattice calculation would also allow for the 
extraction of other phenomenologically interesting quan- 
tities such as magnetic and electric charge radii and mag- 
netic moments. 



2.1 Lattice Techniques 

On the lattice, we determine the form factors Fi{q^) and 
^2(9^) by calculating the following matrix element of the 
electromagnetic current 



2mN 



u{p, s) , (1) 
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where u{p, s) is a Dirac spinor with momentum p and 
spin polarisation s, q ^ p' — p is the momentum transfer 
with — ~q^, TOjv is the nucleon mass and is the 
electromagnetic current. 

The form factors of the proton are obtained by using 



( \ 2_ 1 - 



while for iso-vector (i.e. proton — neutron) form factors 
= "Tm" - dj,,d . (3) 
It is common to rewrite the form factors Fi and F2 as 



q 



(2mAr)2 
F2{q^), 



F2{q\ 



(4) 

(5) 



which are known as the electric and magnetic Sachs form 
factors, respectively. 

At zero momentum transfer, Fi(0) = C/e(0) gives the 
electric charge (e.g. 1 for the proton), while 



1 



(6) 



gives the magnetic moment, where F2^\o) = n^^^ is the 
anomalous magnetic moment. 

In order to extract the non-forward matrix elements 
from our lattice simulations, we compute ratios of three- 
and two-point functions 



7^(^,r;p',p;0) = 



Cr{t,T;p',p, O) 



C2{t,p') 

C2(T,p')C2{t,p')C2{t~T,p) 



(7) 



C2{T,p)C2{t,p)C2{t-T,p') 

which for large time separations, <C r ^ t < \Lt, where 
Lt is the temporal extent of our lattice, is proportional 
to the matrix element we are interested in, {p'\Oq\p). The 
nucleon two- and three-point functions are given, respec- 
tively, by 

C2(t,p) = Tr [ru„poi(S(T,p)B(0,p))] , 
Cr{t,T-p\p,0) ^Tr [r{B{t,p')O{T)B{0,p))] . (8) 

Here t and r are the Euclidean times of the nucleon sink 
and operator insertion, respectively, p ' (p) is the nucleon 
momentum at the sink (source) , and O is the local vector 
current 

0{t) = ^(r)7^V(T) , (9) 

which we renormalise non-perturbatively |10| . The trace in 
Eq. ((HJ is over spinor indices and the F matrix determines 
the polarisation of the nucleon with -Tunpoi = 5(1+74)- 
We note here that in the calculation of nucleon matrix 
elements, we neglect contributions coming from discon- 
nected quark diagrams as these are extremely computa- 
tionally demanding. Hence, in the following we mainly re- 
strict ourselves to the calculation of iso-vector matrix ele- 
ments where the disconnected quark contributions cancel. 

Finally, we use the Sommer parameter, tq, to set the 
scale with tq = 0.5 fm. 
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Fig. 1. \/Q^F2/Fi form factor ratio on three datasets with 
the same pion mass (~ 550 MeV), but with different lattice 
spacings, a = 0.085, 0.080, 0.068 fm. 

2.2 Results 

Of particular interest is the need to understand the be- 
haviour of the form factor F2{Q'^). The question arises 
which is the best way to fit the form factor since such a 
fitting function also allows an extrapolation of the form 
factor to = 0. This is a necessary ingredient to find 
the anomalous magnetic moment of the nucleon, k. 

Based on perturbative QCD, Fi should scale asymp- 
totically as l/Q"^ while F2 - [111112 ■ It is difficult 
to obtain lattice data with high enough precision over a 
large enough range of values to distinguish between a 
dipole or tripole behaviour. It may, however, be instruc- 
tive to consider the form factor ratio F2{Q'^) / Fi{Q'^) since 
asymptotically this ratio should scale as Spin po- 
larisation experiments have instead found that the data is 
compatible with 



1 



F2iQ') 

Fim ™ 



(10) 



To investigate the asymptotic behaviour of the form 
factor ratio F2{Q^) / Fi{Q^) , we plot in Fig. ^ the results 
for \/Q^F2 /Fi obtained at three working points with ap- 
proximately the same pion mass, but with different values 
of the lattice spacing. Here we observe the lattice data to 
be consistent with a constant for > 1.5 GeV^, similar 
to the experimental data. Multiplying these results by an 
extra factor of y/O^, as suggested by perturbative QCD, 
would clearly destroy the plateau. Quantitatively, though, 
the lattice data is higher than the corresponding experi- 
mental ratios, cf ^Jl- This shows that the lattice simula- 
tions are able to reproduce the qualitative features of the 
experimental data, but for a quantitative reproduction the 
pion mass is still unrealistically large. 

In the following we fit Fi and F2 with a dipole ansatz 



Ft\q') 



F^{0) 



(1 - q^Mfy ^^^^ 

where F^^^O) = 1, F2'"^(0) = and M, is the fitted 
dipole mass for the form factor, i. 
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Fig. 2. Results for the isovector magnetic moment as a func- 
tion of m^. The experimental value is denoted by the star. 



We display our results for the isovector magnetic mo- 
ment in Fig.Elas a function of m^. Our results are in good 
agreement with recent quenched |14l[T5ul6, | and Nf = 2 
[TH] results, which indicates that there appears to be lit- 
tle effect due to quenching on the magnetic moments, as 
predicted in (17) . The experimental value is indicated by 
a star at the physical pion mass. We clearly see that a 
linear extrapolation would miss the experimental point. 
This, however, is not completely unexpected as results 
from chiral perturbation theory suggest that we should 
observe a dramatic increase in the results at lighter pion 
masses |14II17| . The new points at lighter pion masses, 

< 0.2 GeV^, are beginning to show a hint of such cur- 
vature, although more work needs to be done to reduce 
the error bars. 

3 Quark momentum fraction, (x) 

Forward matrix elements (no momentum transfer) pro- 
vide moments of quark distributions in some scheme, S, 
at some scale, M: 



(7V(p)|0^^--^">|iV(p))'5 
where 



(12) 



D 



■D 



D = ^{D — D) and {•••} indicates symmetrisation of 
indices and removal of traces. 

Matrix elements with no momentum transfer are deter- 
mined from a simplified version of the ratio of three-point 
to two-point correlation functions given in Eq. See 
[TS] for additional details. 

We use non-perturbative renormalisation as outlined 
in Section 5.2.3 of JS| to convert our lattice results to the 
MS scheme at = 4 GeV^. 

In the language of the parton model, is often de- 
noted by (x"^' 



Jo 



Fig. 3. Isovector (a;) as a function of in the MS scheme at 
fi'^ — 4 GeV^. These preliminary results are obtained at four 
different lattice spacings (in fm): 0.092 (triangles), 0.085 (dia- 
monds), 0.080 (circles) and 0.068 (squares). The star indicates 
the phenomenological result of the MRST analysis j20) as given 
in |18| . This is in agreement with a recent higher order analysis 

m 



Of particular interest is the first (n — 2) moment, 
w| — (x)^, which determines the fraction of the nucleon's 
momentum carried by the quark, q. This quantity is noto- 
rious on the lattice for producing values much larger than 
phenominologically accepted results. These discrepancies 
can possibly be explained by the fact that all lattice calcu- 
lations to date have been performed at quark masses that 
are much larger than the physical masses [19| . Hence, it is 
a challenge for current lattice simulations to calculate (x) 
at small enough quark masses in order to search for the 
severe curvature predicted in Ref. 

Figure 01 displays preliminary results for (a;)'-"^'^^ with 
pion masses as low as ^ 320 MeV. Before we can draw 
any conclusions on the behaviour at small quark masses, 
we need to study scaling violations and finite size effects 
more carefully. Indeed, it has been suggested P^ll^ that 
a volume of at least (4 fm)'^ is required to confirm the 
predicted chiral curvature. 



(13) 4 Generalised parton distributions 



dx I 



{q{x) + {-\Tq{x)\=vl. (14) 



4.1 Matrix Elements And Moments of GPDs 

For a lattice calculation of GPDs, we work in Mellin-space 
to relate matrix elements of local operators to Mellin mo- 
ments of the GPDs. The non- forward matrix elements of 
the twist-2 operator in Eq. (fT!^ specifies the (n — Vf^ 
moments of the spin-averaged generalised parton distri- 
butions. Replacing 7 with 757 leads to moments of the 
spin-dependent GPDs. In particular, for the unpolarised 
GPDs, we have 

dxx''-'H'^{x,tt)^H^{^,t),j' dxx^-'E" ^El, 

(15) 
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Fig. 4. Generalised form factors A"o"'*, A^^"^, A^^"^ together 
with a dipole fit. AU form factors have been normahsed to 
unity. 

where ^ r.-i j 

i=0 

EUtt) = ^n,2.(i)(-2e)'^ - C«(i)(-20"|„ even ■ 

i=0 

(16) 

Here we denote the invariant of the momentum transfer by 
t = = {p' — . The generalised form factors A'^ 2i(*)' 
^n2ii^) ^^"i (^ni^) foi' the lowest three moments are ex- 
tracted from non-forward nucleon matrix elements of the 
operators in Eq. (O ■ 

For the lowest moment, Aiq and Biq are just the Dirac 
and Pauli form factors Fi and F2, respectively 

J ^dxH''{x,tt)^AUt)^F,{t), (17) 

j ^dxE^{x,i,t)^ Bfo it) = F2 {t) , (18) 

while Aio and Biq are the usual axial-vector and pseu- 
doscalar form factors, respectively 

j'^dxH'^{x,tt)^AUt)^9A{t), (19) 

j ^dxE'i{x,^,t)^Bl^{t)^gp{t). (20) 

We also observe that in the forward limit (t = ^ = 0), 
the moments of Hq reduce to the moments of the unpo- 
larised parton distribution ^„o(0) = 

4.2 Results For Generalised Form Factors 

Burkardt 1^ has shown that the spin-independent and 
spin-dependent generalised parton distributions (x, 0, t) 
and H{x^ 0, t) gain a probability interpretation when Four- 
ier transformed to impact parameter space at longitudinal 
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Fig. 5. Generalised form factors A^^''', A^q'^, ^30''' together 
with a dipole fit. All form factors have been normalised to 
unity. 

momentum transfer ^ = 

'^^""^ = / e-^^^-^^Hix, 0, -Al) , (21) 

(and similar for the polarised Aq{x, b±)) where q{x, bj_) is 
the probability density for a quark with longitudinal mo- 
mentum fraction x and at transverse position (or impact 
parameter) b±. 

Burkardt ^ also argued that H{x,0,—Al) becomes 
Z\^-independent as a; — > 1 since, physically, we expect the 
transverse size of the nucleon to decrease as x increases, 
i.e. lim^;^! g(x, b_L) oc 6'^{b±). As a result, we expect the 
slopes of the moments of H{x, 0, —A^) in A^ to decrease 
as we proceed to higher moments. This is also true for the 
polarised moments of H{x, 0, — Al), so from Eq. (|16|l with 
^ = 0, we expect that the slopes of the generalised form 
factors Ano (t) and Ano {t) should decrease with increasing 
n. 

In Figs. 01 and El we show the ^-dependence of Ano{t) 
and Ano{t), respectively, n — 1,2,3, for /3 = 5.40, Kgca — 
Kvai = 0.13500. The form factors have been normalised to 
unity to make a comparison of the slopes easier and we 
fit the form factors with a dipole form as in Eq. (|ll|l . We 
observe here that the form factors for the unpolarised mo- 
ments are well separated and that their slopes do indeed 
decrease with increasing n as predicted. For the polarised 
moments, we observe a similar scenario, however here the 
change in slope between the form factors is not as large. 
The flattening of the GFFs Ano{t) has first been observed 
in Ref. 7 , where at the same time practically no change 
in slope was seen going from A2o{t) to A3o{t). 

Although fitting the form factors with a dipole is purely 
phenomcnological (see Rcf. 24 for an alternative ansatz), 
it does provide us with a useful means to measure the 
change in slope of the form factors by monitoring the ex- 
tracted dipole masses as we proceed to higher moments. 
We have calculated these generalised form factors on a 
subset of our full complement of (/3, k) combinations and 
have extracted the corresponding dipole masses. Recall 
that Aio is the Dirac form factor Fi, while Aiq is the 
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Fig. 6. The lowest three moments of the GPD H{x,(, = 0,t) 
(top) and H{x,(, = 0,t) (bottom) in impact parameter space 
as a function of impact parameter, b. 

axial form factor gA- Hence the dipole fits can be com- 
pared with experiment. A hnear extrapolation produces a 
result larger than experiment for both the polarised and 
unpolarised case, although the findings of Ref. suggest 
that the chiral extrapolation of the dipole masses of the 
electromagnetic form factors may be non-linear. 

In Fig. we show the lowest three moments of the 
GPD H{x,£, = 0,t) (top) and H{x,£, = 0,t) (bottom) 
in impact parameter space. The curves correspond to the 
Fourier-transformation of our dipole ansatz Eq. with 
the dipole masses extrapolated linearly to the chiral limit, 
to 6_L-space, and the shaded error band is a result of the 
errors in the extrapolated dipole masses at the physical 
pion mass. The curves have been normalised so that they 
represent line densities with J dbq"{b) = 1. The top figure 
of Fig. clearly shows how the u — d quark distribution 
narrows as we proceed to higher moments n and thereby 
larger values of the average momentum fraction, while for 
the polarised case in the bottom figure, the narrowing of 
the distribution is not so severe. 
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